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Abstract 

Being based on V. Konoplev's axiomatic approach to continuum mechanics, the paper broadens its 
frontiers in order to bring together continuum mechanics with classical mechanics in a new theory of 
mechanical systems. There are derived motion equations of 'abstract' mechanical systems specified 
for mass-points, multibody systems and continua: Newton-Euler equations, Lagrange equations 
of II kind and Navier-Stokes ones. 

Quasi-linear constitutive equations are introduced in conformity with V. Konoplev's definition of 
stress and strain (rate) matrices. 
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I. INTRODUCTION 



Classical mechanics is based on the axiom system introduced by I. Newton p]. In result of 
generalizations made by L. Euler it is also used to studying the kinematical and dynamical 
behavior of physical objects modeled as a rigid body or their aggregates. 

In the case of a parcel of air, water or rock consisting of a large number of particles, a 
corresponding discrete model, which can be constructed with the help of classical mechanics 
methods, would be hopelessly complicated. A different sort of models has been developed 
over the last three centuries to describe such physical systems. The model, called continuous 
medium or continuum, exploits the fact that in air, water and rock nearby particles behave 
similarly. The corresponding theory discounts the molecular structure of physical systems 
and regards matter as indefinitely divisible (here particles are characterized by their place 
volume and mass density). Thus the intent is to obtain a mathematical description of the 
macroscopic behavior of physical systems rather than to ascertain the ultimate physical basis 
of phenomena. 

The analysis of the behavior of physical systems modeled as a continuum consists that we 
know as continuum mechanics. 

A new architecture of mechanics is suggested in [21 [3j under the conditions that 

1. there are no boxes or particles which can be rotated and deformed; 

2. there are no mass-points (points with zero volume and non-zero mass). 

The first condition makes it essentially various w.r.t. conventional continuum mechanics 
while the second condition deepens the conflict between classical mechanics and that of 
continua [1]: 

. . the dynamics of a continuous system must clearly include as a limiting case 
( corresponding to a medium of density everywhere zero except in one very small 
region) the mechanics of a single material particle. This at once shows that it 
is absolutely necessary that the postulates introduced for the mechanics of a con- 
tinuous system should be brought into harmony with the modifications accepted 
above in the mechanics of the material particle 

Following [21 [3l [5] we aim to remove the conflict by bringing together the continuum me- 
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chanics of Konoplev and the classical one into a theory of mechanical systems. 

The new theory gives mathematical foundations to mechanics, which can be called Newto- 
nian as it remains true to the principles of classical mechanics [H [5] such as the absolute 
space and time, the concept of a mechanical system consisting of points in 3-dimensional 
space as well as those of the mass additivity, actions-at-a-distance and differential laws of 
motion, Galileo's principle of relativity, etc. (it sounds curiously, but I. Newton has defined 
mass, as well as force, as 'the measure of the same' - see, e.g., definitions I, II and VI in [1] 
- as though he has foreseen application of the measure theory to mechanics in 20th century 
013115]). 

The principle demand to a theory of mechanical systems is that 'the problem of mechan- 
ics comes to describing motions being in nature, namely, to their description in the most 
complete and simple form' [B] . Within the framework of this understanding the key concept 
of our theory is that a mechanical system is a set of points equipped with some fields: the 
mass, force, velocity ones, etc. 

The theory is built on relatively simple, transparent ideas, some conventional notions are 
used, but sometimes their sense is radically changed. We try to give all of them on tabula 
rasa without using any background in the field of mechanics. That is why no prior knowledge 
of continuum mechanics or the classical is required. It does not mean that we have done all 
our best in order to avoid any mechanical reminiscences. However giving no comments or 
motivations, we are about to point out all technical details of the introduced constructions 
(for, as Goethe has told, 'God is in the small things ...'). 

To demonstrate the new theory effectiveness we define the main classes of mechanical sys- 
tems and deduce sufficiently many results known in the conventional mechanics: kinematics 
equations and Newton-Euler and Lagrange equations, stress-strain relations, etc. 

We shall use the expression 'see also' in the case where a given statement differs in details 
from that of cited works and thus it is formally absent in them. 

II. PARADIGM OF NEWTONIAN MECHANICS 

The three laws of motion were first compiled by Sir Isaac Newton in his work 'Mathematical 
Principles of Natural Philosophy', first published on July 5, 1687 (in Latin 'Philosophiae 
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Naturalis Principia Mathematica'). Newton used them to explain and investigate the motion 
of many physical objects and systems. For example, in the third volume of the text, Newton 
showed that these laws of motion, combined with his law of universal gravitation, explained 
Kepler's laws of planetary motion. 

Newton's Laws hold only with respect to a certain set of frames of reference called Newtonian 
or inertial reference frames. Any reference frame that is in uniform motion with respect to an 
inertial frame is also an inertial frame, i.e., Galilean invariance or the principle of Newtonian 
relativity. 

Newton's first law is a restatement of the law of inertia. It apparently occurred to several 
different natural philosophers and scientists independently. Aristotle had the view that all 
objects have a natural place in the universe: that heavy objects like rocks wanted to be at 
rest on the Earth and that light objects hke smoke wanted to be at rest in the sky and the 
stars wanted to remain in the heavens. He thought that a body was in its natural state when 
it was at rest, and for the body to move in a straight line at a constant speed an external 
agent was needed to continually propel it, otherwise it would stop moving. 

The 17th century philosopher Rene Descartes also formulated the law, although he did not 
perform any experiments to confirm it. 

The first law (the inertia law), in less accurate form, has published still Galileo, and Newton 
gave credit to him. Galileo, however, realized that a force is necessary to change the velocity 
of a body, i.e., acceleration, but no force is needed to maintain its velocity. This insight leads 
to Newton's First Law - no force means no acceleration, and hence the body will maintain 
its velocity. 

In the given interpretation mass, acceleration, momentum, and (most importantly) force are 
assumed to be externally defined quantities. This is the most common, but not the only 
interpretation: one can consider the laws to be a definition of these quantities. 

Besides, Newton has definitively buried the representation which has taken roots from an- 
tique times that laws of motion of terrestrial and heavenly bodies are absolutely various. 
In its model of the world all Universe is subordinated to the uniform laws supposing the 
mathematical formulation. 

Newton's laws were verified by experiment and observation for over 300 years, and they are 
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excellent approximations at the scales and speeds of everyday life. Newton's laws of motion, 
together with his law of universal gravitation and the mathematical techniques of calculus, 
provided for the first time a unified quantitative explanation for a wide range of physical 

phenomena (here we do not discuss the relativistic mechanics). 

Thus Newton's merit is the decision of two fundamental problems. 

1. Creation of an axiomatic basis for mechanics which has actually passed this science in 
the category of strict mathematical theories. 

2. Creation of dynamics which connects behavior of a body with characteristics of external 
influences on it (forces). 

The Mathematical Principles of Natural Philosophy (Principia) is Newton's fundamental 
work in which it has formulated the law of universal gravitation and Newton's three laws 
being the base of the classical mechanics. 

Newton defined main concepts - mass, force, inertia ('congenital force of a matter'), quantity 
of motion, etc. The absoluteness of space and time which measure does not depend on 
position and speed of the observer are postulated. On the basis of these accurately certain 
concepts (they are a part of Newton's axiom system without which it is wrong.) three laws 
of Newtonian mechanics are formulated. While physicist Aristotle asserted that speed of a 
body depends on motive force, for the first time Newton made the essential amendment: 
not on speed, but on acceleration. 

Newton's original Latin was translated quite closely by Motte (1729). We shall give some 
definitions and the laws with the help of this translation. 

Definition I. The quantity of matter is the measure of the same, arising from its density 
and hulk conjunctly. 

Definition II. The quantity of motion is the measure of the same, arising from the velocity 
and quantity of matter conjunctly. 

Note that the last notion is not quite clearly used at Descartes before. 
Now the quantity of motion is known as an impulse (of a body). 

Definition III. The vis insita, or innate force of matter, is a power of resisting, by which 
every body, as much as in it lies, endeavours to persevere in its present stale, whether it be 
of rest, or of moving uniformly forward in a right line. 
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This force is ever proportional to the body whose force it is ; and differs nothing from the 
inactivity of the mass, but in our manner of conceiving it. A body, from the inactivity of 
matter, is not without difficulty put out of its state of rest or motion. Upon which account, 

this vis insita, may, by a most significant name, he called vis inertia, or force of inactivity. 

Definition IV. An impressed force is an action exerted upon a body, in order to change its 
state, either of rest, or of moving uniformly forward in a right line. 

This force consists in the action only; and remains no longer in the body, when the action 
is over. For a body maintains every new state it acquires, by its vis inertice only. 

The vector nature of impressed forces is discovered in the parallelogram law (see Corollary 
II of the Principia - p. 84 in Motte's translation). 

Definition V. A centripetal force is that by which bodies are drawn or impelled, or any way 
tend, towards a point as to a centre. 

Definition VI. The absolute quantity of a centripetal force is the measure of the same 
proportional to the efficacy of the cause that propagates it from the centre, through the spaces 
round about. 

Definition VII. The accelerative quantity of a centripetal force is the measure, of the same, 
proportional to the velocity which it generates in a given time. 

Definition VIII. The motive quantity of a centripetal force is the measure of the same, 
proportional to the motion which it generates in a given time. 

Law I. Every body perseveres in its state of rest, or of uniform motion in a right line, unless 
it is compelled to change that state by forces impressed thereon. 

Newton's first law postulates presence of such phenomenon, as inertia of bodies. Therefore 
it also is known as the inertia Law. Inertia is the phenomenon of preservation by a body 
of speed of motion (both on size, and in a direction), when on a body no forces or the 
vector sum of all operating forces operate (that is equally effective) is equal to zero. To 
change speed of motion, on a body it is necessary to work with some force. Naturally, the 
result of action of identical forces on size on various bodies will be various. Thus, say that 
bodies possess inertness. Inertness is a property of bodies to resist changing of their current 
condition. The inertness size is characterized by mass of a body. 

It is necessary to notice that Galileo supposed free motion not only on a straight line, but 
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also on a circle (it is visible, from astronomical reasons). Galileo has also formulated the 
major principle of the relativity, which Newton has not included in the axiomatic system 
because for mechanical processes this principle is a direct consequence of the equations of 

dynamics and (sec the Principia): 

CoroUciry V. The motions of bodies included in a given space are the same among them- 
selves, whether that space is at rest, or moves uniformly forwards in a right line without any 
circular motion. 

Newton considered space and time as absolute concepts, uniform for all Universe, and has 
obviously pointed out in the Principia. 

Prom the modern point of view, such formulation of Law I is unsatisfactory. Pirst, 'body' it 
is necessary to replace the term on 'a material point' as the body of the final sizes for lack 
of external forces can make and a rotation motion. Second, and this main thing, Newton 
in the work leant against existence of absolute motionless frame , that is absolute space 
and time (the modern physics rejects this representation). On the other hand, in any (we 
will tell, rotating) frame the inertia law is incorrect. Therefore the Newton's formulation 
requires specifications. 

Law II. The alteration of motion is ever proportional to the motive force impressed; and is 
made in the direction of the right line in which that force is impressed. 

The second law states nothing about the innate force of matter as by definition IV only 
impressed forces are actions exerted upon a body, in order to change its state, either of rest, 
or of moving uniformly forward in a right line. 

It is impossible to consider the first law as a special case of the second one where the motive 
(impressed) force is absent as the former postulates existence of inertial frames while the 
latter is formulated already in such frames. 

The vector nature of the second law addresses the geometrical relationship between the 
direction of the force and the manner in which the object's momentum changes. Before 
Newton, it had typically been assumed that a planet orbiting the sun would need a forward 
force to keep it moving. Newton showed instead that all that was needed was an inward 
attraction from the sun. Even many decades after the publication of the Principia, this 
counterintuitive idea was not universally accepted, and many scientists preferred Descartes' 



7 



theory of vortices. 

In a case when the mass of a material point time-invariant in due course, Newton's second 
law is usually formulated with use of the acceleration notion. 

Newton's second law - the differential law of the motion describing interrelation between 
force applied at a material point and acceleration turning out from it of this point. Actually, 
Newton's second law enters mass as a measure of display of inertness of a body in a chosen 
inertial frame. 

Some authors interpret the first law as defining what an inertial reference frame is; from 
this point of view, the second law only holds when the observation is made from an inertial 
reference frame, and therefore the first law cannot be proved as a special case of the second. 
Other authors do treat the first law as a corollary of the second. The explicit concept of an 
inertial frame of reference was not developed until long after Newton's death. 

Law III. To every action there is always opposed an equal reaction: or the mutual actions 
of two bodies upon each other are always equal, and directed to contrary parts. 

The Third Law means that all forces are interactions between different bodies, and thus that 
there is no such thing as a unidirectional force or a force that acts on only one body. If body 
A exerts a force on body B, body B simultaneously exerts a force of the same magnitude on 
body A - both forces acting along the same hue. Underhne that these forces are enclosed 
to different bodies that is why at all are not compensated. 

Newton's laws, strictly speaking, are fair only in inertial frames. If we fairly write down 
the equation of body motion in a non-inertial frame it will differ from Newton's second law 
by the form. However it is frequent, for consideration simplification, enter certain fictitious 
'force of inertia' and then the motion equations correspond in a kind very similar to the 
Newton's second law. In mathematical relation it is correct, but from the view point of 
mechanics it is impossible to consider the new fictitious force as something real, as result 
of some real interaction. Once again we will underline: 'force of inertia' is only convenient 
convention in order that motion laws appear the same in inertial and non-inertial frames. 

Not all motion equations can be founded in the framework of Newtonian mechanics described 
above. For example, we need the principle of constraint release in the case of non-free 
(constrained) bodies. 
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The sense or senses in which Newton used his terminology, and how he understood the 
second law and intended it to be understood, have been extensively discussed by historians 
of science, along with the relations between Newton's formulation and modern formulations. 
The modern interpretation of Newton's laws is given in [7\. 

With the definitions and laws given above one must connect the definitions and the laws of 
Newton with the concepts of classical mechanics [8] such as the concept of a body (medium) 
consisting of points in 3-dimensional space as well as those of the mass additiveness, action- 
at-a-distance, Galilean (inertial) frame and Galileo's principle of relativity, principle of 
Newton's determinacy, principle of release from constraints and so on. 

It is important to note that all definitions and laws given above are stated without any 
mathematical expressions. Below we shall try to give them with the help of contemporary 
mathematics. 

III. MAIN NOTIONS AND PRINCIPLES OF NEWTONIAN MECHANICS 

In what follows we shall use Galilean spacetime [9] introduced as a quadruple G = 
{A4, V4, g, r} where 

1. V4 is a 4-dimensional vector space, 

2. T : V4 — 7- Ri zs a surjective linear map called the time map, 

3. g = (■, ■) is an inner product on ker{r} (= R3), and 
4- A4 is an affine normed space modeled on V4. 

Introduce a point-wise spatial set X3 with the translation space V3 of 3-dimensional (free) 
vectors and a parameterization t G Ri of the image of r being in a point-wise time set 
T with the translation space V_|_ of 1-dimensional (free) vectors having one and the same 
sense. For some parameterization t G Ri of T the differentiable map Ri — ?■ X3 is called 
motion. 

A. Screw space 

It is considered as conventional fTO] that the screw calculus is not adapted for the description 
of continuous media, and being very attractive representation of a system of forces 
and rigid body motions with the help motors and screws, nevertheless it has no essential 
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practical value . . .' [TT]. As a result in mechanics there is mainly absent the fundamental 
understanding (concept) that the interaction between mechanical systems is described with 
the help of screws. 

The using of the screw concept is the key for the theory of mechanical systems (including, 
a continuum, a mass point and a rigid body) which is below constructed. 

Let us define two fields and Tl of vectors attached to points of X3 such, that for any two 
points a and b G X3 with r^, r;, G 9^ and fla, fib & ^ there is the following relation 



def 



ra = n = r, Ha = fib + ab X r 
where x means vector product. 



(1) 



Definition 1. [12j The set {fH, 9H} is called screw while {f,fla} is element of reduction of 
the screw at a point a G X3. The vectors f and fla are called main vector (resultant) and 
total moment of the screw (at the point a), respectively (the screw moment is not that of the 
vector f as we do not connect the point a with a point in X3/ 

We do not support the idea to use the name 'torser' from the French word 'torseur' instead 
of 'screw' 1121. 



A screw with the property f x fla = (for all points a G X3) is called slider [T2j. We say 
that a slider is axial at some point a, if /2a = (we have at least one such point). It is useful 
to note that the sum of axial sliders is the axial slider, too. 

A screw with the property r = is called couple, from ([l]) follows that fla = flb'= fl for all 
points a and b G X3 [T2] . 

Definition 2. We shall call wrench and twist the following operator forms, respectively: 

„wr 

defined at the point a G X3. 

We may also define them as objects with the following properties: 
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where / and O are unit and zero matrices. 

Henceforth we shall briefly say slider a wrench with the slider property and note it as V . 

Define the triple of orthogonal unit vectors e° = {C]^", Cg", 63*^} in the 3-dimensional space 
V3. Let us introduce 6 wrenches (twists) such that at the point a G X3 their elements are 
defined as follows 




where o G V3 is the null vector. 

As any screw is defined in the unique way by its element of reduction at some point, these 
6 wrenches (twists) generate the basis c** of the screw space, the first triple of the wrenches 
being axial sliders and the second one being couples. As elements of the spaces V3 and R3 
are called vectors, we may use the names wrench, twist and slider for coordinate columns of 
elements of any wrench, twist or slider in the basis e" at the point a G X3. 

A screw can be resolved in a sum of a slider and a couple if it is neither slider nor couple 
[12] . This resolution is not unique. A couple can be represented as sum of two sliders. That 
is why any screw (as a vector in the screw space S) is a slider sum, too. This fact is used 
[2] in order to define a slider as the primary notion of screw theory. 

Remark 1. For a given system of line vectors there exists a point a G X3 such that its main 
vector 7" and total moment Ttj = are such that 1^ x 74 = 0. In this case the system is 
called screw being the set of the following elements [TS^ 

— the straight line (screw axis) passing through the point a; 

— the main vector 7" giving the screw axis sense; 

— the moment vector Tti (being collinear to 

One says that the point a is that of the screw reduction and this screw depicts a screw motion. 

This definition leads to no matrix tools of screw calculus which may simplify the reduction 
of line vectors to the simplest equivalent system fT^ . 
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B. Main measures of Newtonian mechanics (see also [2113]) 

Let us define the Lebesgue measure fit on a-algebra of subsets in T while on a-algebra of 
subsets in X3 there be so called determinative time-invariant measure 

where /iac(^) is the absolutely continuous component w.r.t. Lebesgue measure /is and fipp{A) 
is the pure point (discrete) component presented as fipp{A) = J2k f^pp(-^k) for points in an 
arbitrary subset A E (73 such that fipp{xk) 7^ 0. These points are called pure, the others 
being called continuous [14j. We assume fiac to be Lebesgue measure /is. 

We shall further use the measure /x^^ for definition of points with mass, but without volume, 
and bodies with volumes, but without masses and forces exerting on them. 

Definition 3. Let A G a^, then the measure m{A) : A ^ Hi is called mass (measure of 
inertia) . 

Due to the Radon-Nikodym theorem [T3] we may specify m{A) as Lebesgue-Stieltjes integral 

"^(^) = J XAPxf^Lsidx) 

with a /i^^-integrable (mass) density px w.r.t. the measure ix^g{dx) (here Xa the charac- 
teristic function of A). The density can be time-varying. 

The set X3 C X3 is called set of concentration of the measure m on X3 if m{B) = for 
everyone /i^,, —measured set i? C X3 \ X3. 

We shall use the notion of signed measure [E] being a generalization of the concept of 
measure by allowing it to have negative values. Some authors call it charge, by analogy with 
electric charge, which is a familiar distribution that takes on positive and negative values. 

Let ri{-) G V3 be a function on whose components (in some basis) are signed measures. 
Then: 

1. the function ri{-) is called vector signed measure on as; 

2. a function C(") ") ^ "^3; defined on x as and being a vector signed measure by each 
of arguments, is called vector signed bi-measure; 

3. the vector signed bi-measure ({■,■) is called skew if ({A,B) = —({B,A) for any A 
and B G a^. 
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Definition 4. Given A and B G o"3, the skew vector signed hi-measure B): {A, S) — )■ S 
is called measure of action of B on A. 

Remark 2. It is important to point out that C. Truesdell defines mainly the measure of 
action as ?>— dimensional vector /J/. 

We specify F{A, B) as Lebesgue-Stieltjes integral 

wliere elements of /x^g-integrable slider functions l^^^'^^ and l^^y'^^ are axial at x G A and 
y & B, respectively (they can be represented as corresponding Lebesgue-Stieltjes integrals 
with densities being axial sliders). 

The set A'^ = X3 \ A is called environment of A. It is clear that 

T{A, A^) = T{A, A^ + A) I = J (2) 
We shall assume that I't'^^^^''^ = on the set X3 \ X3. 

Definition 5. The slider function /'^(^'^'') is called intensity of the action of upon x E X3. 
Remark 3. The intensity can also depend on the motion prehistory. 

Exemplify the introduced notion. Let the skew bi-measure g{A, B) G V3 be such that 




where 7 is a positive (gravitational) constant, elements of the /i^g-integrable slider function 
l^i^'^") are axial at x G X3. 

Definition 6. [3] The slider function pj^^^'^"^ is called intensity of gravitating action o/x^ 
upon X G X3. 

C. Fundamental principles of dynamics 

Let at be Borel cr-algebra of subsets in T while as is Borel a-algebra of subsets in X3. 

Let us fix some parameterization t G Ri of T, then the differentiable bijection: X3 — Xt C 
R3 is called motion, t is a time instant. 
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For any point x G X3 the motion defines the point x{t) G Xf. Introduce the radius-vector 
f^it) = (Oo, x{t)] called position of x{t) G and the vector = r^it)' called its velocity 
w.r.t. Oq. Thus we equip the set X3 with the fields of positions, velocities and the measures 
of mechanics. 

Let the slider l^"^ be axial at x G X3. 

Second Newton's law (see also [31 [16]). There exist a Cartesian frame So with the origin 
Oq and a parameterization t G Ri of T such that motion of a point x & A'^ is described by 
the following equations in the slider form 

1. if the point x is continuous 

p.{i:-'y = e-^^^'" (3) 

2. if the point x is pure 

mM:-'y = f^pp{x)it^'''">'' (4) 

where = pxf^pp{x) is mass of the pure point (coordinate representations of vectors 
in So are marked with the superscript ° while the subscript means that the slider 
moment is computed w.r.t. the point x; to honor Newton, we use the superscript " for 
derivatives by t). 

Henceforth we call the parameterization and the frame So Galilean (this formulation of 
second Newton's law is connected with first one and isolated systems nohow). 

Remark 4. In the case of time-varying densities of inertia (masses) relations (^^^ are 
invariant w.r.t. Galilean group while the traditional form of second Newton's law 

does not. In such case relations (^^@ include in themselves slider functions of so called 
reactivity (see the well-known equation of Mescerskii). 

Definition 7. The set G fK3,T), (0-3, at) and {fij^g,m, J^) answering the second Newton's 
law and the principles of causality, determinacy and relativity /5/ are called Universe of 
Newtonian mechanics, elements of a— algebra cr^ being called mechanical systems f^. 

In the given definition (see also [21 E]), similarly to that of probability space [17J, Universes 
of mechanics are separately specified for every mechanical problem. 

From relations (|3|-(|4]) follows that motion of a point x G A'^ is described in the following 
slider form 
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1. if the point x is continuous 

p.{ir'y = it'""^^'' (5) 

2. if the point x is pure 

mAir'y = f^ppi^)iP'''^^^'' (6) 

It is easy to see that all slider forms of motion lead to the following vector relations 

1. if the point x is continuous 

p^v"- = /ipp(x)0°(x, x'') (7) 

2. if the point x is pure 

mxV°- = yUpp(x)0°(x,x^) (8) 

In the motion equations the intensities are defined nohow, and any action intensity pictures 
some mechanical system [9] in depending on its 'constitution'. Sometimes some part of the 
intensities is implicitly given, while another one must be defined from the restriction or con- 
straint imposed on a point, its velocity and, perhaps, derivative of the velocity, beforehand 
set, i.e., not dependent on the law of point motion. In this case a point which motion is in 
agreement with constraints is called constrained. 

Example. Let the vector x^) describe the action of x'^ on x when constraints are absent 
and the constraints be given by the equation (T(r°, t>°, t) = where a is a differentiable 
vector-function of the instant t, the position r°. and the velocity v^. After differentiating 
(T(r°, t>°, t) = we have 

(9cr n da ri da 
t'. + -^v^- + — = 



arS'^ " ■ dt 

and 

,,0. = ( )Ti 9a , da ^ , da , da 

OVx OVx OVx OTx Ot 

if the above inverse exists. 

Hence there exists such slider function l'^'^^^^"^ (axial at x G A'^) that 

p^(/--0)- = /J(---^)'0 + /f'"'^'° (9) 

if the point x is continuous or 

mx{Q''y = /xpp(x)[/j(^'^^)'° + /f •^^)'°] (10) 
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if the point x is pure. 

In this way we may introduce the following 

Principle of constraint release. Motion of any constrained point x & A'^ is described by 
equations ^ and (10) (in the Galilean frame Eq) with some jji^g-integrahle slider function 



j^c{x,x'') qqIIq^ intensity of constraint action upon x G A^. 

The principle of constraint release demarcates two categories of actions, namely, active and 
passive ones: it says that an active (motive) action creates motion while a passive one only 
puts obstacles in this motion. If we remove constraints then only active actions are kept. 

Remark 5. One must not suppose that the principle of release from constraints and that of 
DAlembert eliminate the difference in the nature of active forces and passive ones (constraint 
actions and forces of inertia). It is only for the sake of convenience that we use these 
principles: only forces the resultant of which is f exert on a point (body) 



IV. A MASS POINT 

Consider a set A G consisting of a unique pure point of the measure /i^g as a free mass- 
point. In this case from relation ^ follows the well-known second Newton's law (in the 
Galilean frame So): 

m.vl- = fl (11) 
where f^ = fipp{x)(j){x,x'^) is called force exerting on the mass-point x. 

Note that in the case where the mass is time-varying the force includes in itself that of 
reactivity (see the well-known equation of Mescerskii). 

Motion of a mass-point can be constrained. Let us give the description of constrained 
motion. 

Variety of constraints contains so called ideal and non-ideal ones. Ideal constraints generate 
constraint actions having the direction and sense of the normal to the corresponding man- 
ifold. We shall assume that constraints are ideal {Axiom of ideal constraints), scleronomic 
and holonomic. 

Ideal holonomic and scleronomic constraints force the point under consideration to move 
along with a certain manifold cr(r°) = having lower dimension than its configuration 
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space. Let this manifold can be parameterized with some vector q. The vector q is called 
generalized one, its first derivative being called generalized velocity q\ 

For any point x of the manifold we have 

= v{<i) 

If the columns of the following matrix [19] 

drj 

are linearly independent, they form a basis of the linear space T = T(g) being tangent to 
the manifold at a point q. 

If the columns of the following matrix [19] 

- 1^ 

are linearly independent, they form a basis of the linear space N = N(g) being orthogonal 
to the manifold at a point q. 

Thus R3 = T X N with the basis [r, u] . It is easy to see that 

P,R3 = T, P.Rs = N 
where Pr = {t'^t)~^t'^ and P^, = v{v^v)''^v^ are projections. 

It is obvious that r°" = rg" and r^" = rq" + {-^q')q' (here the matrix ^^q' is square). 



As the constraint is supposed to be ideal and therefore PtT^ = 0, from (10) follows 

dr 

^.[q" + Pr{Q^q-)q-]=Prfl (12) 



Applying the projection Py to (10) we define the following equation 



(9r 

m^Py{-^q-)q- = P.(/° + c^) = P./° + c° 
From this relation follows that 

dr 

c° = m^Py{—q-)q- - P^f^ 

i.e., the constraint force is not a function of time, but it depends on the generalizing coor- 
dinates and velocities as well as on the active force. 

It is clear that the theory above can be applied to mass-point systems. 
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V. RIGID BODIES 



Definition 8. (see also [20]). A bounded closed set A & cr^ is called rigid body if 

1. constraints applied on its points keep distances between them not changing with time; 

2. the constraints are ideal. 

A rigid body may contain continuous and pure points. 

Remark 6. In elementary manuals of mechanics, transition from a mass point to a body as 
a point system is made somehow imperceptibly; constraint forces are not mentioned at all, 
and instead of a lawful exception there is an illegal, silent exclusion of these forces. They 
remain ordinarily without any attention and even without a mention, as if they did not exist 
at all. 

Multiplicative groups of motions. In V3 let us have 3 triples of orthogonal unit vectors 
e° = {ci", 62", 63"}, e*' = {el ,62,6^} and e'^ = {6^,62 , e^}, where p and k are naturals 
(here we may use p = 0, p = k — 1 and p = k). 

Assume that e*^ is chosen as the main basis and the space R3 has so called canonical basis 
e°, consisting of = col{l,0,0} and 6° = col{0,l,0}, 6° = col{0,0, 1}. The other bases 
can be movable. With their help let us define the frames Sq = ^(Oq, e°), Sp = S{Op, e^) and 
Sk = £{0k, e'^) with the origins Oq, Op and Ok- 

Define rotation matrices Co,p and Cp^k such that Co^pCp^k = C'o.fc and for any free vector A 
there are the following relations 

A*^ = Co,pA^, A^ = Cp^k^^ 
where A*^, A^ and A^ G R3 are the coordinate columns of the vector A in the bases e°, e^ 
and e*^. Hence we have also A° = Co,fcA'^. 

Let X be an arbitrary point fixed in £p. Introduce the radius-vectors tq and Tp G V3 of the 
point x w.r.t. the origins Oq and Op, respectively. Define (io,p = tq — fp. Then we may 
represent the relation tq = do^p + rp G V3 in £q as Tq = d^p + Co,pr^. As is time-constant, 
with differentiating the last relation we have f ° = v^ p + C^ pr^ where f ° = Tq and Vq p = d^'p 
are velocities of x and Op w.r.t. Oq in the frame £^05 respectively. Hence 
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For any vector / 



G R3 introduce the cross product matrix 



Let us define (in Sp) [2T] : 



r 



dcf 



-fs f2 
-/2 /l 



(13) 



the coordinate column d^^ G R3 of the translation vector of £p w.r.t. £q] 

the coordinate column t>Q G R3 being known as quasi-velocity oi the translation 



of Sp w.r.t. £q] 



— the cross product matrix Wp^p =^ Cp_oC'6,p where the triple cJq j, G R3 is known as 
angular quasi-velocity of rotation of Sp w.r.t. £0 and is the eigenvector of Co,p, 
answered with the eigenvalue 1; 

— quasi-velocity V^^p = co\{vq p,ijjQp} G Re of motion of Sp w.r.t. Sq. 

These algebraic quantities are answered with geometrical ones, e.g., vectors of the translation 

velocity vo,p G V3 and the instantaneous angular velocity cJo.p G V3 which are defined with 

the help of the basis e^. We may use them in order to define the twist (kinematical screw) 

V with the element Vo,p = coljwo.p, i^o,p} at the point Op. 

In the kinematics the angular velocity ujo,p defines the rotation axis of Sp. 

Introduce following matrices 



rpO 

0,p 



MX 



o 
I 



0,p 



O Cq^p 



(14) 



wr,0 



L^^pTi^'^'P where vr^'^ and 71^^ are 



Theorem 1. [3] Let us have a wrench vr""". Then tt\ 
its elements at the points Oq and Op, the matrix Lq ^ has the representation 



T wr 
^0,p 



^0,p-^0,p 



(15) 



and belongs to the multiplicative group C^^'(Jl,6) such that 



T wr. T wr ^wr /hwr 

^0,p — ^0,p^O,pi ^0,p 



px 
0,p 



O 



px px 
^0,p <p 



(16) 
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Proof. The representation of Lq^ follows directly from the screw definition. 



Relation ^ is true as from 1^ follows that L^;- = Tq^Cq^p + T^^pC^^p = T^,pC^,p{C^,oC^;, 

i/^® T-iO. N T wr (f,wr 

"'"'-"p.O-^O.p'-'O.p/' ~ -'^O.p^O.p- 

The matrices of the kind Lq ^ form a group because there are Lq^L^j^ = T^ pC^^pT^ i^Cpf, = 

7^0 7^0 r^® _ 7^0 _ Twr c cnhinHpY n anrl T^^-^ — fT^ r"® "l^l — C'^''^ (T^ 

^o,p'''p,k^Q,p^p,k — ^o,k^o,k — ^o,k ^ suomaex p ana i^q ^, — l,Jo,fc'-"o,feJ ~ '-'o,A: 1-^o,A:J 



~ '-"fc,0-'fc,0'-"0,fc ^k,0 — '''kfi^kfl ~ ^k,0- 

The similar statement is true for twists vr*'": tTq""'" 



Lq'^tt*"''^ where we have the matrix 



T tw 

^0,p 



T tw ffitw fhtw 
-^O.p^O.p' ^0,p 






/ 







I 







T wr 
^0,p 







/ 


I 



belongs to the multiplicative group C*'^{Tl, 6) such that Lq"^' 



wr,T 
0,p • 



Note that in contrast to the groups of motions in the 3-dimensional space the groups 
£""'(7^,6) and C*^{TZ,Q) are multiplicative. 

Newton— Euler equation. Let the frame £p be attached to a body Ap. 
Lemma 1. [3J There is the following relation 



]Vx ,p 

'p 



Proof. The statement is true as 



'p 



^p ^0,p^ ^p 



p 



px 



^px\2 



I 




I 




yP = 




^px 




^px 


_ p 




_ P 



H,p + ^o.p^p) 



I — r 



r-PX 
P 

px\2 




where the relation Wq.p'^p ~ ~'"p^'^o,p used. 



Due to the lemma we have (C^)" = {L^^pi;-^)- = L^:p{e^Vi;p+ ^Z^^V^p) or 



;^p^o!p + ^o!p^p^o!p)- That is why from (|9|)-(|l0| follows that 



where (and henceforth) all integrals are taken w.r.t. Lebesgue-Stieltjes measure /i^^; the 
set Ap is immobile in the frame £p. 

According the rigid body definition the constraints are considered as ideal and thus 

/x.,^P^o/f'^^^'V.sW = 
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Theorem 2. [3] The motion of Ap (w.r.t. £q in the frame Sp) is described by the (Newton- 
Enler) equation (see also fl^ ) 



^0,p ^ ^0,p ^p ^0,p ~ 



(17) 



where Op = JxAp^pPxl-i'Lsi^^)> ^ ~ IXaJo^^'^ ^'^f^Lsi^^) '^^ wrench calculated in Sp 
and generated by the main vector and the total moment acting on the body Ap. 

Systems of consecutively connected bodies [22]. Let us consider a system of A; + 1 
consecutively connected bodies Ap, p = 0,k (the body Aq is immobile). Its motion is 
depicted by the following Newton-Euler equation 



AV- + BVa = Fa 



(18) 



where A and B are known matrices, Va = coI{Vqp}, Fa = co1{Fq }, p = l,k. 
Newton-Euler equation (18) is considered w.r.t. 'absolute' quasi-velocities p of the bodies 



(calculated in Sp w.r.t. the main frame Sq). But in practice there are only the 'relative' 
quasi-velocities of the frame Sp w.r.t. Sp-i. Thus we must connect the 'absolute' 

quasi-velocities with 'relative' ones. 



Lemma 2. [3J For a system of consecutively connected bodies there is the following compo- 
sition rule 

J^O,p J \ "^O.p-l / \<^p-l,p 

Proof. As the bodies are connected consecutively there is the relation Co,p = Co,p-iCp_i,p is 



Cp^p-iUQ J_i+ul_i . Besides de 



true. With differentiating it we have loq.^ 

fine the vectors cTo.p-i = (Oq, Op-i] and rfp_i,p = (Op_i, Op), then d^p = c?o,p_i + rfp_i,p, v^p 



^0,p^ 



^1 + '^p-i,p5 '^p-i,p ~ C*o,p-ic?p_i p, dp'_ip — Vp^ip + C'Qp_idp_ip — Co^p^iUiQp_idp_ip + 



-'p-l.p "p-l,p'^0,p~l 



>-l,p ' "p,p-l'^0,p-l- 



I P— 1 X ^ ^P— 1 

'p-l,p — ""p-l.p + '-^0,p-ll^o,p-l'-"0,p-l"p_l 

nence u^p — f^o,p-i ^ ^p-i,p ^ ";p,p-i^p-i,p5 ^o,p ~ "-"p.p-i ^o,p-i ^ '-^p,p-i"p,p-i^o,p-i ^ 

Cp^p- 1 o 



^0,p 

O Cp^p-i 
Hence we have 



-' "p,p-i 



O 



YP-'^ I T/P _ TtW T/P-1 I T/P 



k=p 



-trP — T tw -trk \/k 

0,p ~ 2-^-'^p,k^k-l,k-> ^k-l,k 



^k-l,k 



(19) 



k=l 



fc-l,fc, 
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(19) 



Va = LVr 



(20) 



where K = col{ V}^, . . . , V^., 



rk-l 



,Vn-i n}y ^ is the trian- 



gular matrix with blocks ^ being functions of 'relative' frame rotations and translations. 
Thus we have 

ALV; + {AL- + B)Vr = (21) 



where L' is analytically calculated due to relation (16). 



It is easy to see that the matrices of relation (21) depend on rotation matrices (and linear 



and angular quasi- velocities, too) that is why equation (21) must be considered along with 
the Euler kinematical relation 



a 



p-i,p 



L'p-i,pi^p-i, 



(22) 





Fig. 1. Multibody system graphs 



Multibody systems with tree— like structure. Consider a multibody system with tree- 
like structure given by the graph in Fig. lA. Let vertices f represent the system bodies or 
the origins of the attached Cartesian frames where the index % numbers the tree-tops, the 
index j numbers the bodies from the base to the corresponding tree-tops. Introduce V^-^ as 
quasi-velocities characterizing rotation and translation of the frames w.r.t. Then we 
have the sets Vl^^ , Vl^, V^t VH}, V,'l V^, V^i, V,]^, V^i}, {V,'[\ V^i, V^i, 
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^3,4 , , ^5,6 }> ^1,2 , ^,3 , ^3,4 I, IK),! ^1,2 ^2,3 ^3,4 }, IK),! . ^1,2 , ^2,3 , ^3,4 I and {V^q,! > 

T.1,1 -r^l,! ylM ryO,l yl,l yl,l yl,l yl,2 y2,2. ryO,l yl,l yl,l yl,l yl,2 y2,3. ryO,l 

K)!2^K)!3^K)!4^}'{^o°l^^o^2^K)!3^''^o^4^^ with the same subscripts as in 

the case of consecutively connected bodies for the relative and absolute quasi-velocities. 



This case is considered above that is why we arrive at relation ( 20 ) with the known matrix 
L and 

ITTaO,! ta1:1 Ta1:1 T^l.l T^l.l Ta1.2 ta2,2 ta2,3 t^1,4 Trl,5 Tr5,5 ta5,6t 



irT/0,1 taI.I Trl,l Trl.l Trl.l 1^1,2 t/^.^ TALI'S TA^j^ taI.O t aOiO TaO,D-| 
Vr = Co\{Vo [ , V;^2 , ^2,3 . ^3,4 , ^4,5 . ^4,5 . ^5,6 , ^5,6 . ^3,4 , ^2,3 . ^3,4 , ^3,4 I 

Remark 7. T/ie results obtained can be immediately applied to systems with loops, e.g., if 
in the system under consideration (see Fig. IB) the vertex 6^ is connected with 6^ by the 



2,2 



r2,3 Trl,4 Trl,5 T^5,5 T/5,6l 



edge (Q ,Q ). In this case relation (17) is the same, but in the case where constraints are 



considered there are the following additional constraints (5^, 6^) + (6^, 6'^) + (6^, 5^) = and 



^2,2^2,3^2,2 T 

'-"5.6'-"6,6'-"6,5 ~ ■ 



Parameterization of rotation matrices. The order of system (21)-(22) may be reduced. 



To this end one uses different parameterizations of rotation matrices. 



Euler angles. Let Cp^i^p = C1C2C3 where 




1 
COS ip —sirup 
sin ip cos ip 



C2 



cos Ip — sin ^|J 
sin cos Ip 
1 



(23) 



cos sin 
1 
- sin -i? cos ip 

are so called the simplest rotation matrices; ip, i), and ip are Euler angles [23] . 

Introduce the triple Ap_i^p = col{ip, -t^yip} as a parameter. Then there is the matrix -Dp-i,p 
such that [2] 



Hence equation (21) must be considered along with the following relation 

\- ^ n-i , ,p 

p—i,p p— i,p p— i,p 



(24) 



(25) 



and Cp_i_p = Cp_i^p(Ap_i^p) if the matrix Dp\^p exists. 

Fedorov vector-parameter. To parameterize rotation matrices we may introduce Fedorov 
vector-parameter 
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Definition 9. [2^ The number triple f G R3 is called Fedorov vector-parameter of a 
rotation matrix C, if it is answered with the following matrix (see ^) 

The inverse map of Cayley restores the rotation matrix 

C=(/ + /X)(/-/^)-l 

It is easy to be verified (for example, by means of Maple®) that the following relations are 
true 

_ (l-||/f)/ + 2//^ + 2/x 



Let the rotation matrices Co,p and Cp^k have Fedorov vector-parameters /o,p and fp^k- It is 
known that they are eigenvectors of these matrices, i.e., 

Co,pfo,p = fo,p, Cp^kfp,k = fp,k £ R-3 

As the space R3 has 3 bases e°, and e'^ we may write /o,p = /o°p = /o,p, fp,k = fp^k = fp,k 
and define the following vector 



fp,i 



i i i 

where col{/i, Z^, /s} = and col{/?, /"} = Co,p/;,,. 

Definition 10. The vector fp^k is called vector o/Rodrigues (the half of the vector of finite 
rotation - see 121^). 

The definition is motivated by the fact that so defined vector is coUinear with the 
instantaneous angular velocity Up^k, and thus it defines the rotation axis. 

Remark 8. In [24J it is explicitly pointed out that vector-parameters are Gibbs vectors which 
can be defined in the form Qp ^ = '^^ficf' (in mechanics free vectors and their coordinate 
columns in the canonical basis e° are known with the same name vector as elements of the 
vector spaces V3 and R3, and there is no other basis except the canonical one in l2^ ). In 
general the vector (jp^k does not coincide with fp^k (see Fig. 2), as the bases e°, e^ and are 
different. 
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Fig. 2. Gibbs and Rodrigues vectors. 
There is the following relation 



As the vectors Up-i^p and rp-i,p are collinear, we have relation (24) where Xp-i,p = fp-i,p 

2 



and -Dp-i,p - i+y/^.i^pP (-^ + /p-i,p) 



Thus equation (21) must be considered along with the following relation 



f;-i,p = ^ (1 + ii/p-i,pir)(/ - fp-i,pni + fp-i,pH-i,p (28) 

and Cp-i^p = Cp_i^p(/p_i^p). 

Euler-Rodrigues parameters. To parameterize rotation matrices we may use quaternions. 
Definition 11. The set A = {Aq G R, A G V3} is called quaternion. 

Quaternions generate the algebra with the quaternion product 

A o M = { Ao/io - (A, p) , Ao/i + /ioA + A x /I} 
where M = {yUo,/^}, is the inner product. 

Any vector A can be imaged as a quaternion A with the zero scalar part. That is why we 
may define the quaternion product of two vectors A and jl as follows 

A o /i = {— (A, /I), A X /I} 
There exists the unit quaternion Ap_i^p = {Aq, \p-i^p} (with ||Ap_i^p|| = 1) such that 

^p-i,p°^p-i,p° ^p-i,p = ^p-i,p^ ^p-i,p = ^2Ap_i_poAp_-|^ p, = 2 '^p-i,p°^p-i,p (29) 
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where Ap_i p = {Aq, —^p-i,p} is conjugation of Ap„i_p. 

Let us denote col{a;i, U2, W3} == ^l-i,p, col{Ai, A2, A3} =^ ^l-i,p and col{Ao, ^l-i,p} '= ^l-i,p 
then the orthogonal matrix Cp_i^p corresponding to a rotation by the unit quaternion Ap_i p 
is given in the following form 



Cp-i,p(Ap_i_p) 



Aq + A^ — A2 — A| 2A1A2 — 2A0A3 



2A1A3 + 2A0A2 



2A1A2 + 2A0A3 Aq — A^ + A2 — A3 2A2A3 — 2AqAi 
2A1A3 — 2AqA2 2A2A3 + 2A0A1 Aq — A^ — ^2 "I" ^3 



(30) 



The quadruple Ap_^ is known as that of Euler-Rodrigues parameters. 



From (29) follows 



—LUi —U2 —UJ3 

OJl UJ3 —UJ2 

U2 —0J3 ui 

U3 UJ2 —OJl 



(31) 



Hence equation (21) must be considered along with relations (|29|)-(|31|) 
From (30) follows also that there is the matrix -Dp-i,p 



i5p_i,p(Ap_]^ ) such that relation 



p-i,p ~ -^p-i,pAp_ip is true. 



Lagrange equation of II kind. Let Ap_i p be a triple of Euler angles or Fedorov vector- 
parameter. 

Definition 12. 1. The vectors gp_i,p = col{dp_i p, Ap_i,p} and q^^i^p = co\{dp'_^ p, -^p-i,p} 
are called canonical generalized coordinates and velocities of the frame £p in the motion 
w.r.t. the frame £p-i; 
2. the relation 



p-i,p 



diag{/, -Dp-i,p} 



(32) 



is called equation of kinematics of £p-frame w.r.t. £j 



p-i- 



From relations (21) and (24) follows the Lagrange equation of II kind 



A{q)q" + I3{q,q-)q- =T 



(33) 



where A{q) = L'^M^ALM, B{q,q-) = L'^M^[ALM- + {AL + S)M], J" = L'^M^Fa, M 
diag{Mp_i,p}, q = col{gp_i,p}. 
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In the many cases there are constraints on motion of multibody systems, and the matrix 
exists such that the matrix N'^N is non-degenerate and we may introduce the generahzed 
coordinate qc = Nq G Rm where the natural number m is not more 6/c [12]. Then from 



relation (33) follows 

Acq'c + -Beg; = J^c 

where Ac, Be and J-'c are known matrices and column. 

As to the quadruple Ap_j^p, we may replace \p-i,p with A^^^^ in the above definition and 



equation (33). It is clear that the corresponding matrix A proves to be singular. Under some 
assumption this equation is equivalent to a system of differential equations in Cauchy form 
and algebraic ones. The algebraic equations can be treated as constraints on the multibody 
system motion. It means that we may introduce 'new' coordinates, e.g., Euler angles or 
Fedorov vector-parameter, in order to obtain the Lagrange equation with a non-singular 
symmetric matrix A. 

Scholium . The singularity of Lagrange description in different generalized coordinates and 
velocities is the price that we must pay if we give up Newton-Euler description in kinematical 
twists. In practice this price is not very high. 



VI. A CONTINUUM 

Suppose the set A G X3 has no pure point of the measure /x^g, ii^g{dx) = ixz{dx) and 
m{dx) = PxfJ^sidx) in A'^. 

Strain matrix and its rate. Given x{t) and y{t) G A'^ in the instant t G T, define their 
radius-vectors rx{t) and ry(t) (in £^0) and the vector h{t) = Vy — Txif). If h{t) is small we 
have 

vl{t)^vl{t) + dvlld{rl)h%t) 
Define the matrix Zx{t) as the solution of the following equation 

Z-At) = dvl/d{rl) 

with initial data = I for t = Iq. 
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Definition 13. [3] The matrices and Z'^ are called strain one and its rate at the point 
x{t) G A'^ in the instant t, respectively. 

There is no reason to consider tlie strain matrix and its rate as important (kinematical) 
characteristics of continuum motion. 

Stress matrix. Let us define (see also [21 E]): 

1. a section S between the set A and an arbitrary plane P; 

2. the vector bi-measure 

V{A) = I x,/^(^'^^V3(rfa;) 

where elements of the slider function are axial at x E A; 

3. on the set 5* the slider function l^^^^ of the measure V w.r.t. Lebesgue 2— dimensional 
measure fi2{dx) on Borel a-algebra a2 of open subsets of S such that 

V{S) = I Xs^'^"'"^V2(rfx) 

4. 3 X 3-matrix-function T^,. of x and t which can be differentiable by x the necessary 
number of times and such that the vector x*^) has the coordinate representation 

6''{x,x^)=T,nl 

where Hx is the normal to the plane P at the point x; 

5. the entries of A°(a;,a;'^) being connected with the rows (j = 1,3) of the matrix 
by the following relation (in the frame Sq) 

A°(x, x') = DivT^ = col{divr^\ divTj, divT^^} 

Remark 9. One may see that the measure T>{A) is introduced under the influence of Gauss- 
Ostrogradsky divergence theorem f^, but here it is said nothing about the properties of T^, 
e.g., about its symmetry. 

Definition 14. [3J 1. The slider function l^^^'^"^ is called intensity of stress action upon 
X G A"; 

2. Tx is called stress matrix. 
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Notion of continuum. A matrix-function of entries of some matrices is called isotropic if 
it is invariant w.r.t. iSO(R, 3). Assume that the matrices T^, and Z'^ are invariant w.r.t. 
rotations. An isotropic map T^lZ^, Z'^) is called constitutive or stress-strain relation. 

Let us note the set of all isotropic maps from the strain matrix Z^ and its rate Z'^ to the 
stress matrix as a'^{Zx, Z'^). 

Definition 15. Suppose that 

j^<p{x,x<') _ pj^g{x,x<') j^c{x,x^) _ ^A(x,x<=) ^g^^ 

the stress matrix belongs to a'^{Zx,Z'^) and the measure of inertia is time-constant on 
A'^, i.e., ^m{dx) = 0. Then the set A is called continuous medium or continuum. 

Motion of continuum. Due to relations ^ and (34) the equation of continuum motion 
at a point x ^ A^ is of the form (see also [31 (in the Galilean frame Sq) 

pxvl = Px9° + Div T^, p'^ + divp^ = 0, = T^{Z^, Z;) e a^{Z,^, Z;) (35) 

where Div Tr is the constraint action |25] . 



Some constitutive relations. Constitutive relations define the properties of a continuum 
and its motion equation. 

3— dimensional case . For any 3 x 3-matrix A the aggregate PAQ is an isotropic function 
of A if the matrices P and Q are proportional to / with scalar coefficients being invariant 
w.r.t. rotations. 

Define the following linear combination 

T = nEi + ra^s + rg^s (36) 

where rj are invariant w.r.t. rotations (they can be functions of the time, invariants of 
and so on), Ux = Z^ or f/^ = Z'^; 

E, = (tr[/,)/, E2 = symU, = 0.5([/, + t/J), E, = antU, = 0.5([/, - Uj) (37) 

Theorem 3. [26] All isotropic 3 x 3-matrix functions of entries of Ux are given by relation 
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Thus due to the theorem we may define the following relation 

= -roEo + ri^i + ra^a + rg^a, Eq = I (38) 
as the most general linear constitutive one. It is conventional the invariant w.r.t. rotations 



Tj to be called rheological coefficients (w.r.t. the set of E.i). The constitutive relation (38) 
is called quasi-linear if its rheological coefficients are functions of [^-matrix invariants (in 
particular it means that the summand riEi could be omitted in quasi-linear constitutive 
relations) . 

If Ux = Zx and tq = the continuum is called elastic material, if Ux = Z'^ and tq > (called 
Pascal pressure) the continuum is called viscous fluid |2Zj • 



Remark 10. Continua defined by relations (31)-(3S) coincide with the continua used in 
continuum mechanics in the following cases [2Z] 

1. the Pascal pressure tq is positive and ri = r2 = = (ideal fluid); 

2. To is non-negative and = (continua of Navier-Stokes-Lame type); 

3. To is non-negative and ritr/ + r2 = — )■ trT^,. = — rotr/ (continua used in some 
theories). 



In order to use relations (37)-(38) in the motion equation (35) we must calculate 

DivTx = — gradro + trf/j.gradri + riDiv (trf/^./) + 0.5r2(Div Ux + Div Uj) + 

0.5r3(Div Ux - Div Uj) + 0.5( Ux + [/J)gradr2 + 0.5( Ux - t/J)gradr3 (39) 

Introduce the next notations 



Ux 



ul 


ul 


ul 


ul 


ul 


ul 


u\ 


ul 


ul 



u 



ik 



du\ 



r^ = col(r?,r0,r3°) 



Let u^ be the rows of Ux, then with the help of routine calculations we have Div(tr?7^/) 
co\{u\^,ul2,ul3} and 



DivUx 



Ull + ^22 + "^33 



ufi + M22 



U 



33 , 



^u\, + ul^ + ul^ 



Divf/, 



T 



u\i + U22 + ^23 
"^31 + "^32 + "^33 
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2— dimensional case . In the case of 2 x 2-matrices it is easy to see that for matrices P and 
Q the aggregate PU^Q is an isotropic map of Ux if P and Q are of the kind al + al where 
-0-1 

/ = , the scalar coefficients a and a are invariant w.r.t. rotations. 

1 

Introduce the following matrices E2 = Ux, -E3 = Ux^ , E4 = lU^, E^ = lUx^, Eq = U^I, 
E-j = Ux^I, Eg = lUxI, and Eg = lUj^I- The linear combinations of the matrices Ei 
{i = 2, 9) generate a manifold with a basis consisting of the linear independent matrices E2, 
Es, E4 and E^. 

We may use this basis in order to define the following relation 

Tx = -ro/ + fo J + ntr UI + npf UI + rs + r^Uj + rju, + rjuj, pfU = tr{7 [/} (40) 

as constitutive relation with rheological coefficients and (w.r.t. the set of Ei) being 
invariant w.r.t. rotations. 

Correct continua . A continuum is called correct if the corresponding constitutive relation is 



invertible |3] (see map (38) and (40)). 

Let us stop at 2— and 3— dimensional constitutive relations with the same 'structure'. With 
the help of routine calculations we see the following statement to be true. 



Theorem 4. In 2- and 3-dimensional cases let constitutive relations be of the ( 38)-form and 
(ritr/ + r2)r2r^ ^ where I is used as 2- and ?> -dimensional identity matrices, respectively. 
Then there exists the inverse map 

Ux = UqI + ni(trr^)/ + n2symr3; + n^aiAT^ 

where 

ro —ri 1 1 

no = , ni = r, n2 = —, n^ = — 

ritrJ + r2 r2(ritr7 + r2) r2 

Thus Navier-Stokes-Lame continua are incorrect as r^ = 0. 

Rheological coefficients and moduli (ratio). The three coefficients 

1 r2(ritr/ + r2) 1 r2 rii ri 



77-2 — ni ri(tr/ — 1) + r2 ' ^ 2n2 2 ' ni — 77-2 ri(tr/ — 1) + r2 
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can be called Young modulus, e, shear or rigidity one, fi, and Poisson ratio, v, respectively 
(see also [3]). Note that there is the known relation e = 2/i(l + u). 

In this way of definition Young modulus, shear or rigidity one and Poisson ratio depend on 
continuum dimensions. 

Let us define now 

Eo = 1, Er = (trf/,)/, E2 = U„ E, = (41) 



then we may take the constitutive relation (38) with new rheological coefficients tq, ri, r2 



and rs (w.r.t. the new set ([4l])) and new Young modulus, shear or rigidity one and Poisson 
ratio. Thus they are not unique. 



VII. BRIEF COMMENTS 

Continuum mechanics is closely connected with Riemann integral theory. In continuum 
mechanics as well as in Riemann theory there is realized the idea of approximating an area 
by summing rectangular strips (segments, squares or boxes), then using some kind of limit 
process to obtain the exact area required. It is safe to say that we may name the well known 
mechanics of continua as that of Cauchy (due to the man who created it). 

The Riemann integral, natural though it is, has been superseded by the Lebesgue or 
Lebesgue-Stieltjes integral and other more recent theories of integration. In this way, V. 
Konoplev suggested a new architecture of continuum mechanics based on Lebesgue integral 
and his algebraic theory of screws. As result in Konoplev mechanics there do not arise boxes 
or particles which can be rotated by the laws of Newtonian mechanics as well as there are 
no imaged surfaces with stresses over them and other concepts of Cauchy mechanics. But 
with introducing the measures as Lebesgue integrals he was forced to exclude mass-points 
and their systems from consideration. 

Unlike V. Konoplev, we use Lebesgue-Stieltjes integral in order to introduce main mechan- 
ics measures and classes of mechanical systems such that mass-points, rigid bodies and 
continua (under the special assumption about interaction in mechanical systems and their 
'constitution'). In this way we become closer to mechanics of C. Truesdell. 
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CONCLUSION 



It is a first attempt to represent elements of Konoplev's axiomatics and its (possibly debat- 
able) modification in the form of a journal paper. One must realize the difficulties and gaps 
issued from this goal. 

It is impossible to separate the theory given above from that of Konoplev. That is why the 
paper author prefers to yield the palm to Prof. V. Konoplev but carries full responsibility for 
all lacks of this paper. This is the place to express his sincere thanks to Prof. V. Konoplev 
for the collaboration of many years. 
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Appendix: A review 

Dear Professor Cheremensky, 

I have looked at your paper, and I find it nearly incomprehensible. I regret to say that I 
cannot consider it for the Archive of Rational Mechanics, as it does not meet basic standards 
of clarity that would allow me to review it. I should add that, even if it could be rewritten 
in a comprehensible way, the subject matter is probably not suitable for the Archive today. 
During a brief period in the 1970s, the Archive did become a forum for some axiomatic 
work in mechanics, but in hindsight this has become some of the least influential work that 
has been published in the Archive, and we do not encourage it now. Generally, we also 
encourage authors to use the most conventional notation and to avoid abstraction for the 
sake of abstraction, so that the work is readable by the widest possible scientific audience. 
Thank you for your interest in the Archive and I hope you find a receptive audience for your 
work. 

Sincerely yours, 

Richard James, Editor in Chief 

The works of C. Truesdell (Founder and Editor of Archive for Rational Mechanics and 
Analysis: 1952-1989), W. Noll and B.D. Coleman, etc., are published in this journal. We 
are not assured that they are really least influential. They are contributed to foundational 
rational mechanics, whose aim is to construct a mathematical model for treating (continuous) 
mechanical phenomena. 

We would be highly grateful with whoever would bring any element likely to be able to 
make progress the development, and thus the comprehension, of the paper. Any comments, 
reviews, critiques, or objections are kindly invited to be sent to the author by e-mail. 
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